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1. Introduction 

Recently, the four-point amplitude in ABJM theory [1] has been computed at two 
loops [2, 3, 4] that is the first non-trivial order where perturbative corrections appear. 
This amplitude shares many remarkable properties with its analogue in M = 4 SYM, 
namely it is dual conformally invariant [5, 6], it exhibits WL/amplitude duality [7, 8] 
and can be consistently thought of as the first term in an exponential resummation of 
the perturbative series analogous to the BDS ansatz in four dimensions [9, 10]. 

Actually, the amplitude itself divided by its tree level expression strikingly looks 
very similar to the four-gluon amplitude in Af = 4 SYM at one-loop divided by its tree 
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level counterpart. Precisely, the expression of the former ratio, evaluated in dimensional 
regularization d = 3 — 2e, exactly matches the latter in d = 4 — 4e, up to a constant 
and after identifying the renormalization scales. 

In Af = 4 SYM theory, the IR divergent part of the amplitude is completely fixed 
by an evolution equation which constrains its dependence on the renormalization mass 
scale to be proportional to the cusp anomalous dimension / SYM (A) [11]. The finite part 
is equally fixed by off-shell dual conformal invariance [6, 8] and the result is an iterative 
expression [10] 



log ^ = [IR div] + log 2 (I) + const. + 0(e) (1.1) 



In particular, when restricted at first order in the A-coupling, this identity gives the 
one-loop ratio on the LHS in terms of the first order expansion of the cusp anomalous 
dimension, as required by the Ward identities at this order. 

Given the matching between the two-loop ABJM amplitude and the one-loop 
M = 4 SYM amplitude, the Ward identities satisfied by the three dimensional am- 
plitude are the same as those in four dimensions. Therefore, for the two-loop ratio 
of the three dimensional theory we expect an expansion similar to (1.1). In fact, the 
explicit result found in [2, 3, 4] can be factorized as in (1.1) and, quite remarkably, 
the coefficient in front of the finite remainder agrees with the ABJM cusp anomalous 
dimension determined through integrability arguments [13] 



/<,(A) = k YM (A) (1.2) 
2 



being h(X) the ABJM interpolating function [13, 14]. 
Therefore, at the order we are working we can write 



A {2) 

^3d 



A tree 



A {1) 

•^Ad 



A tree 
, ^4d 

JCS 



+ const. + 0(e) (1.3) 

/SYM 

where the mapping (1.2) between the two cusp anomalous dimensions is meant. 

Assuming that this identity can be uplifted to all orders, we are tempted to con- 
jecture that both amplitudes exhibit an iterative structure related by 



log 



A tree 
"^3d 



, Am 
= log 



A tree 
, ^Ad 

JCS 



+ const. + 0(e) (1.4) 

/SYM 



Since the identification is up to (9(e) terms, the perturbative series for the two ampli- 
tudes will not coincide. In fact, when exponentiating eq. (1.4) to obtain the complete 
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amplitudes the (9(e) terms on the RHS will mix with the e-poles, spoiling the identifica- 
tion of the amplitudes order by order. The exact mapping between the two amplitudes 
can be reconstructed only once we know all the 0(e) terms. It is then important to 
check the validity of eq. (1.4) to all orders in e. 

In Af = 4 SYM it happens that in the four-point amplitude the divergent and the 
non-constant finite parts in (1.1) are completely captured by the one loop contribution 
in a very precise way, encoded in the well-known BDS ansatz [10]. Assuming that also in 
ABJM theory the first non-vanishing contribution to the four-point amplitude dictates 
all the non-trivial dependence on the kinematic invariants through its e expansion, it 
is sufficient to find an all-order-in-e relation of the type (1.4) only for the lowest 
non-trivial order in A corrections to both amplitudes. 

In this paper we carry out such a program by proving analytically an exact identity 
between the two-loop ABJM amplitude and the one-loop Af = 4 SYM amplitude to all 
orders in e. This is achieved by first acting on the four dimensional amplitude with a 
differential operator. The result can be manipulated in two different ways, so deriving 
two differential equations. In the first one, the RHS features the 3c? amplitude. In the 
second one, a six dimensional integral appears which can be re-expressed in terms of 
the original 4d amplitude, by means of an integration by parts identity and Passarino- 
Veltman reduction. Equating the RHS of these equations yields a close, all-order-in-e 
relation between the ABJM and Af = 4 SYM lowest order amplitudes. Up to order e 2 
we have managed to check explicitly this identity by evaluating the three dimensional 
momentum integrals up to 0(e 2 ) and comparing the result with the known expansion 
of the 4d integrals at that order [10]. 

The powerful identity we have found allows to rewrite the BDS-like ansatz for the 
ABJM four-point amplitude in terms of the four dimensional one. It follows that, if 
the conjectured ansatz is correct, the similarity between the two amplitudes uncovered 
at first order will propagate all over the perturbative expansion. 

Assuming the exponentiation ansatz to hold in three dimensions, we can speculate 
on the form of the four-loop amplitude. The essential ingredients determining the 
non-trivial parts of the four-loop correction are the cusp anomalous dimension which 
is known up to four loops [15], and the expansion of the two-loop amplitude up to order 
e 2 . In particular, the cusp anomalous dimension along with the two-loop amplitude 
at order zero in e, allows to fix the finite remainder of the amplitude, as shown in 
[3, 4]. In this paper, exploiting the knowledge of the subleading terms in e, we make 
an almost complete prediction for the amplitude at four loops, up to constants and 
scheme-dependent coefficients appearing in front of simple poles. 

The plan of the paper is as follows. In Section 2 we review the result for the ABJM 
four-point amplitude at two loops and discuss its similarity with the four dimensional 
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one-loop amplitude. In Section 3 we first work out the expression of the two-loop 3d 
amplitude up to order e 2 , by explicitly solving the corresponding momentum integrals. 
This allows to realize that the matching with its four-dimensional cousin persists up 
to order e 2 . Motivated by the observation that there must be at least a technical 
explanation for this surprising similarity we then give the explicit derivation of an all- 
order-in-e identity between the two quantities (see eq. (3.27)). In Section 4, assuming 
the validity of a BDS-like ansatz for the 3d amplitude, we propose an all-loop relation 
between the two amplitudes. Finally, in Section 5, assuming the exponentiation ansatz 
to be valid, we work out an almost complete prediction for the ABJM four-point 
amplitude at four loops. A final discussion and five Appendices with all technical 
details follow. 



2. ABJM four— point amplitude: A review 

In the ABJM theory, the four-point amplitude for two scalars and two fermions has 
been calculated at two loops in [2, 3, 4] . Its explicit expression divided by the tree 
level counterpart reads 



(2e)* (2 £ ) 



,,r wr* (^r! + i l0g2( f )+ifi+o(e) 

2.1) 

where A = N/K is the ABJM 't Hooft coupling (see Appendix A for conventions), // 
is the IR scale of dimensional regularization, conveniently redefined such as to absorb 
the e -1 pole 

// 2 = 87re- 7E n 2 (2.2) 

and K\ = 4^ + 3 log 2 2 is a numerical constant. 

Working in M = 2 superspace and using an ordinary diagrammatic approach, the 
result (2.1) arises by summing contributions from six super-Feynman diagrams, after 
performing D-algebra reduction and computing the corresponding momentum integrals 
in dimensional regularization, d = 3 — 2e, up to order 0(e). We list the results for the 
relevant integrals, referring to [4] for a detailed explanation of their origin. 

Keeping the notation close to the one used in Refs. [3, 4], the two-loop amplitude 
can be written as 

Mf d = (Ait) 2 [J(°> (s) + 1^ (s) + 6/ (d) (s) - 2I& (s, t) + (s «-> *)] (2.3) 

where 
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Integral (a) 



[ , ~ (4ir)iP(l-2 £ ) UJ 1 ' 



Integral (6) 



_ o r(i/2 + e )r 2 (i/2-6)r(i + 26)r 2 (-2e) 

1 {S) ~ 2 (47r)T(l-2e)r(l/2-3e) [~7 J { } 



Integral (d) 



Iid)(s) = r3(i/2- e )r(i + 26)r^(- 2e ) 

i; (47r) d r 2 (l-2e)r(l/2-3e) V s / 



Integral (/) 



/</> ( S t) = (l + s/t)T%l/2-e) 

1 ' ' (4vr) d r 2 (l-2e)r(l/2-3e)(t/ / u 2 ) 2e 1 ' ' 



+ioo 



x / ^.r(-v)r(-2e-v)r*(-i-2 e -v)r 2 (i + v)r(2 + 2 e + v) f-V 

J 2iri \t / 



The last expression is the result of using the Mellin-Barnes (MB) representation of the 
Feynman parametrized integral. The star stems for shifting the integration contour 
on the right of the first pole of T(— 1 — 2e — v), ensuring a well-defined expression. 
This integral can be evaluated in the e — > limit [3, 4], leading to the finite result 
(|log 2 S /t + 3C 2 ). 

The expression (2.1) is nicely akin to the one-loop four-point amplitude in M = 4 
SYM theory [16, 17] 



A S YM^V(e) = A 

SYM 



(^r (t/^r , i ^fs\ + 2^_ 



e 2 e 2 



+ iv(f) 



+ 0(e) (2. 



arising from the evaluation of a single box integral. Here, the 't Hooft coupling is 

8tt 2 

// 2 = Ane-^ n 2 (2.9) 



defined as A SYM = irr an d the regularization scale is 



Neglecting terms that vanish when removing the IR regulator, the two expressions (2.1) 
and (2.8) are indeed identical up to numerical constants, once we identify the scaling 
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parameters and shift e — > 2e in M\ d loop {e) to take into account the different loop order. 
At this stage we can then write 

M { 2(e, //) = M^(2e : //') + const. + 0(e) (2.10) 

A stringent question which arises is whether the identification between the two results 
holds at any order in e. The main motivation for investigating this problem comes 
from the observation that in M = 4 SYM all-order terms in the IR regulator are 
crucial for determining the correct exponential resummation of scattering amplitudes 
[10]. Therefore, an answer to this question may shed some light on the structure of the 
exponentiation of scattering amplitudes in the ABJM model, as we now explain. 

Compelling evidence suggests that the four dimensional result (2.8) is the first 
order expansion of an exponential resummation of the perturbative series [10] 



M id = exp 



EM /(0(e) ^ (k)+c(0(e; 

.1=1 



(2.11) 



where Mf] is the four-point one-loop amplitude to all orders in e divided by the 
corresponding tree level expression, while C^(e) contain constants plus (9(e) terms. 
The expansion is in powers of the dimensionless effective coupling, while the mass scale 
(2.9) is hidden inside the one-loop amplitude. The functions /^(e) have an expansion 
in e 

f il \t) = if + f?t + f?t 2 (2.12) 

whose zero order terms coincide with one-quarter the coefficients appearing in the 
perturbative expansion of the scaling function / SYM arising in the dispersion relations 
for magnons 1 . 

In the ABJM theory, perturbative contributions to scattering amplitudes can occur 
only at even powers of A = N/K due to the invariance of the theory under the discrete 
symmetry K — > —K, V ■<->■ V, A ■<->■ B. Therefore, the two-loop contribution (2.1) is 
the first non-trivial quantum correction. 

This observation accompanied by the impressive similarity (2.10) hints that the 
all-order amplitude for ABJM might equally exponentiate as in the four dimensional 
case. In other words, the amplitude (2.1) could be the first order expansion of 



Mm = exp 



oo 

J> 2 < (F l He)Mfhle)+CW(e)) 
.1=1 



(2.13) 



1 The scaling function is found to be twice the cusp anomalous dimension T cusp which controls the 
UV divergences of Wilson loops near the cusps. 
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where the J\f = 4 SYM coefficients /^(e) have been replaced by their ABJM counter- 
parts, /« = jf + ~f?e + /fe 2 . 

A first suggestive support to this ansatz comes from the observation that, according 
to our two-loop calculation, at this order the coefficient = matches the ABJM 
scaling function obtained through integrability in a rather independent context 

fW = q/sym(A S ym)I v^^ >h{x) (2.14) 

with h(X) being the interpolating function appearing in the dispersion relations for 
ABJM magnons. 

Further support to the ansatz (2.13) should come from higher order-in- A results 

(2) 

for which the knowledge of the A4 3 J amplitude at subleading orders in e becomes 

(2) 

mandatory. Moreover, according to eq. (2.10), at two loops we can trade M 3 J(e) in 
eq. (2.13) with M^(2e). If the identity (2.10) were to persist at higher orders in e 
we could express the BDS-like ansatz for three dimensional amplitudes in terms of the 
four dimensional one. Therefore, it is important to evaluate the amplitude at finite e 
and investigate whether the identification (2.10) holds at any order. 

3. An identity between ABJM and M = 4 SYM four-point am- 
plitudes 

In this Section we discuss the evaluation of subleading-in-e contributions to the four- 
point amplitude in ABJM theory. 

By direct inspection of 0(e 2 ) terms we find a refinement of eq. (2.10) which holds 
at that order. Thereafter, we prove a general identity between the ABJM and the 
M = 4 SYM amplitudes at all orders in the regularization parameter. 

3.1 0(e 2 ) identity for the two-loop amplitudes 

When evaluating the subleading-in-e terms of the two-loop ABJM amplitude, it is 
more convenient to work with the IR scale fi" in eq. (2.9) rather than // in eq. (2.2). 
In fact, this allows to avoid the appearance of non-trivial functions of the kinematic 
invariants multiplying powers of log 2 which would make the comparison with the four- 
dimensional amplitude more obscure. The price that we pay is the emergence of a simple 
pole divergence which modifies the relation (2.10) as 

Mf](e, //') = M^(2e, //) - {s^ + t" 2e ) 1 log 2 + const. + 0(c) (3.1) 
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To check the consistency of this relation beyond O(e ) terms, we write 



(3.2) 



j=-oo 



and determine the coefficients A_ 2 , ■ ■ ■ ,A 2 by explicitly computing the integrals (2.4- 
2.7) up to e 2 . The result is reported in eq. (C.2). 

There, we also give the one-loop amplitude of Af = 4 SYM, eqs. (C.3, C.4), where 
we have chosen to write the result in the same form as eq. (3.2) and a doubling of the 
customary regularization parameter e has been performed (d = 4 — 4e). 

Quite remarkably, the expressions for the A_ 2 , ■ ■ ■ ,A 2 coefficients of our amplitude 
are almost carbon copies of the C- 2 , ■ ■ ■ ,c 2 coefficients of the one-loop amplitude in 
M = 4 SYM. A closer look reveals that most of the differences are due to additive 
numerical constants which depend on the subtraction scheme that we choose. 

The only non-trivial difference between the two sets of coefficients is the appearance 
of a log 2 (s/t) term at order e 2 . However, it is easy to see that at least at this order we 
can absorb it in a scheme redefinition 



where D(e) is given in (D.4). This result shows that the connection between the 3d 
and the 4d amplitudes persists at order e 2 . 

In the next Section we prove that this connection is not accidental but can be 
extended to all orders as an exact identity that we derive analytically. 

3.2 All order identity for the two— loop amplitudes: An analytical derivation 

First of all, we observe that writing the amplitude as in (3.2) and taking into account 
the particular dependence of the integrals (2.4-2.7) on the Mandelstam variables, the 
only non-trivial contribution to the coefficients Aj comes from the integral In 
fact, all the other integrals produce just constant factors. 

Therefore, for the time being we concentrate on In the following, we are going 
to prove an exact relationship between this integral and the four dimensional amplitude 




(3.3) 



thus obtaining the following empirical relation 



Mf} (e, M = M$ (2e, //) + D(e) + 0(e 3 ) 



(3.4) 




(3.5) 



where 




(3.6) 



8 



and B(e) is a constant given by 



r*(-2e)r(l + 26) 27Be 



r3(i-4 e) ^ < 37 > 

At second order, fi^ coincides with the expression found in (3.3). 

To prove the identity (3.5) we start with the all-order-in-e expression of the integral 
as given in eq. (2.7). Neglecting the mass scale (we will recover it at the very end 

of the derivation) and defining x = s/t it reads 

, _ W*r'(l/2-e) 
^ 7rr(l/2-3e)Hl-2e) t 1 + ' 

— r(-w)r 2 (-2e - w)r 2 (i + v )r(i + 2e + «) x v 

(3.8) 

At the same time, we consider the Mellin-Barnes representation of the <id amplitude, 
again neglecting the mass scale. Being it given by a single box integral I^\s,t), we 
can write 

M$(2e) = -8n 2 stl£\s,t) 

2ft~ 2e r f dv 

' -T(-v)T 2 (-l-2e-v)T 2 (l+v)T(2 + 2e + v) x v 



(4tt)- 2 T(1 -4e) J 2ni 

(3.9) 

The contour of integration in this expression is ill-defined in the e — > limit, signaling 
the emergence of e-poles. Therefore, we extract the divergent contributions by suitably 
deforming the contour. This leads to 

M$ = r(1 ^ 2 ' ) _H7 2 ' ) f 26 (log(x) + 2^°>(-2e) + 7e - ^ (0) (1 + 2c)) (3.10) 



2(47r)- 2 T(-4e) 



^(i^Firr^) / i^Vi - 2e - + ,)r(2 + 2e + „) 



where ip(°\x) is the digamma function defined in (B.8) and the Mellin-Barnes integral 

that survives has a well-defined contour in the e — > limit and contributes to the 

amplitude beginning at order e. 

t (D 



The product M.^ s 2e depends only on the ratio x = s/t. Deriving with respect to 



x we obtain 



d_ (1) _ r(i + 2e)r 2 (- 2e ) i 

^^4^ J" 2(47T)-2 T (_ 4e ) X 

- 2er(1 _ 4e) y ^-^^ - ^ + + 2e + ^ 



(47r)- 2 T( 
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where the same Mellin-Barnes integral as in (3.8) appears. 
Thus, comparing (3.11) to (3.8) we can write 

(1 + x)^(M^s^) = -2c (An)H [T(s) + T(t) t*] + j^f e M{ d (3.12) 

where A(e) is given in eq. (3.6) and we have defined 

T(s) = TV + 1 m) = r(l + 2e)T\-2e) 1 

{ ) (4tt) 2 - 2 T(1 - 4e) s 1+2e ' K > (4vr) 2 - 2 T(l - 4e) t 1+2e 1 J 

Note that T(s) and T(t) are one-mass triangle integrals in four dimensions in the s 
and t-channel, respectively. 

Equation (3.12) establishes a differential relation between the Ad amplitude and 
the contribution to the 3d one. In order to obtain an algebraic relation, we derive 
a first order differential equation relating the four dimensional box diagram to itself. 
This can be done by using an algorithm similar to the one of Ref. [18]. 

We consider the Feynman-parametrized form of the Ad box integral in d = A — 4e 
introduced in (3.9) 



^(l,!,!,!;*,*) = 



d 4 ~ 4e k 1 



(2vr) 4 - 4e k 2 (k - Pl ) 2 (k - Pl - p 2 ) 2 (k + p 4 ) 2 

1 W- - . (3-14) 



(4^)2-2, t 2~2e J L J (a ia3a; + a2a4 )2+2e 



where the measure is [da] = da\da2da^ S(^2i oti —1). The labels in 1^ indicate the 
powers of the propagators according to the order in which they appear. 
By taking the derivative with respect to the ratio x = s/t we obtain 

xi^t/f) = (1 + 2e)s^tlP - r(3 + ^f + " / [da] 



dx y 47 v ' 4 (47r) 2 - 2e J 1 ( ai a 3 x + a 2 a 4 ) 3+2e 

(3.15) 

The second piece of this equation is proportional to a six-dimensional box integral with 
two indices raised by one unit. Precisely, 

i 



(4vr; 



3-2e 



/ ^ „ TY^-e = 4 1} (2,l,2,l; g ,t) (l + 3 + )4 1} ( g ,t) (3.16) 

J (ai«3 s + a 2 tt4rj + 



where n 4 " are the operators which raise and lower the power of the n-th propagator by 
one unit. 
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Comparing (3.16) with (3.15) we obtain 

x-^(s 1+2e 1 1^(3, t)) = (1 + 2e)s 1+2£ t I ( l\s, t) - Air s 2+2 H (1+3+) /^(s, t) (3.17) 

We further manipulate (1 + 3 + )Iq~\s, t) with the scope of re-expressing it in terms 
of four dimensional integrals. Applying Integration-by-parts relations arising from the 
identity 

n = I d6 ~ 4ek d ( (k- Pl - P2 r \ 

J (4ir) e -^dk» \k 2 (k - Pl ) 2 [(k - Pl - p 2 ) 2 } 2 (k + p 4 ) 2 J V ' } 

we are led to 

4tt s (l + 3+)4 1) (s, t) = 4tt(1+ + 2+ + 3+ + 4+)/ ( f ) (s, t) + Atx Ae 3+J^ (s, t) (3.19) 

Using the Feynman-parametrized form of the first four terms, it is not difficult to 
ascertain that by the condition ^\ = 1 imposed by the delta-function we obtain 

4tt(1+ + 2+ + 3+ + 4 + )4 1) (s, t) = (s, t) (3.20) 

For the last term, comparing its Feynman-parametrized form 

^'W) = *$±£ I w (aia3S ;; a4f)2+ , (3.2D 



with the Feynman-parametrized form of a four-dimensional vector-like box integral 

d 4 ~ 4e k fc" 



h ~i (27T 



)4-4, fc2(A . _ pi )2 (A . _ pi _ p2 )2 (A; + p4 ) 5 



T(2 + 2e) /^ Q ] Q 2Pi + a 3 (Pi + P2Y - o^j ^ ^ 



(4vr) 2 - 2e J L J (aia 3 s + a 2 a 4 £) 2+2e 


we find that it coincides with the (pi + p 2 )-direction of the vector-like box integral 
in four dimensions. This component can be easily evaluated by employing Passarino- 
Veltman reduction and we obtain 

An Ae 3+/« (s, t) = ^-/f ( S , t) + Ae (3.23) 

where T(s) and T(t) have been defined in (3.13). 

Collecting the results (3.19, 3.20, 3.23) and inserting back into eq. (3.17) we obtain 
the desired differential equation for the four-dimensional box integral. Recasting it in 
terms of the amplitude M^(2e) = Sir 2 stl^\s,t) we finally have 

(1 + X) Tx {M ^ {2e) S2€) = + 2e ^ft^int) - T(s)} (3.24) 
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Comparison with the RHS of (3.12) produces an algebraic equation relating M^ d to 



M, 



(i) 



id 



M( d (e) = A(e) (M«(2e) + (4tt) 2 (sT(s) + tT(t))) 



(3.25) 



Finally, reinserting the scale parameters and absorbing the A(e)(47re _7E ) 2e factor as in 
(3.6), this identity casts into the form (3.5). 

In order to recover the whole 3d amplitude in terms of we add to (3.5) the 
contributions from the integrals I^ a \ 1^ and 1^ where we apply the /i^-scheme re- 
definition. Given that 



(47T) 2 



A(e) (47re-TE- 



2c 



(&\s) + I^(s) + 6I (d \s) + (s <* t)) 



= -( s - 2e + r 2e ) 



e 2 ^ e r(l-6e)r 2 (| -e)r 3 (l-e) 
47rr 3 (l - 2e)r(l - 3e)r(l - 4e)r (| - 3e) 



r 2 (-2e)r(2e + 1) 



erf \ -e] -2r(i-2e)r Q + e 



+r 



i 



3e T 2 



1 



e r 2 



+ e 



= ( S - 2e + r 2e ) E(e) 



we finally obtain 



(3.26) 



(i) 



(3.27) 



where we have defined D(e) = (s~ 2e + t~ 2e )(B(e) + -E(e)). Series expansions for A(e), 
5(e), D(e) and 5(e) can be found in Appendix D. It is easy to check that neglecting 
subleading terms for e — > we are back to the relation (3.1). 

4. Three dimensional BDS ansatz revisited 



The identity (3.27) enables us to reformulate the BDS-like conjecture for the ABJM 
four-point amplitude in terms of the original M = 4 SYM all-loop proposal. 

In fact, reformulating the ansatz (2.13) in the fi A scheme, first of all we can write 



Mz d {e,ii A ) = exp 



oo 

£A 2 < (fW(e)MV{k,» A ) + cW(ej) 



(4.1) 
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where the functions fcs\^) are the /^(e) counterparts in the /i^-scheme. Their 
leading coefficients f^ are still determined by one-quarter the Chern-Simons scaling 
function, as the change of scheme affects only f^ l \ 

Modifications in the constant part of the amplitude due to the scheme change are 
included in the new coefficients Ccs\e). 

The convenient choice of the /z^-scheme allows to use the identity (3.27) in the 
previous expression, thus leading to a suggestive ansatz for the all-loop four-point 
amplitude in 3d 

M M {e,n A ) = exp 

= M M 

where we have defined 

n«\e) = fjg>(e) D(le) + C^(e) - C«(2e) (4.3) 

and C^(2e) are the functions appearing in the 4c? BDS ansatz (2.11). 

It is important to note that, apart from a factor (s~ 2U + t~ 2le ) hidden inside D(le), 
the coefficients T-L^(e) contain only constant terms, while the hard-core of the ampli- 
tude is completely encoded in M.^- 

The ansatz (4.2) reveals a deep intertwining between the conjectured exponentia- 
tion of four-point amplitudes in M = 4 SYM and ABJM theories. Therefore, we are 
led to conjecture that the remarkable connection uncovered at lowest order [2, 3, 4] will 
propagate over their entire perturbative series. 



oo 

£ (/£°00 M$(2le, //) + fg\e) D(le) + CgP(c)) 
.1=1 



l/(e)->/cs(e) 



exp 



X> 2 ^«(e) 



l=i 



(4.2) 



5. A conjecture for the four— loop amplitude 



In spite of this beautiful result, the poor knowledge of the functions / C s( e ) spoils the 
power of (4.2) in predicting higher order corrections to the four-point function. Never- 
theless, using known results for the scaling function of the ABJM theory, we are able 
to formulate an almost complete prediction for this amplitude at four loops in terms of 
the M = 4 SYM amplitude at two loops. 

The main ingredients for carrying out this program are the expression 



fg>(e) (M^(2e^") + D(e) 
+fif(e) (-Mi 1 j(4e,/i") + D(2e) + cfi>(e)) 



(5.1) 



obtained from the expansion of the ansatz (4.2), and the value of the ABJM scaling 
function at four-loops [15] 



/ CS (A) =4A 2 -24C 2 A 4 + G(A 6 ) (5.2) 

which leads to 

fcs = fcsfl = 1 /cs|o = — 6C2 (5.3) 

Unfortunately, while integrability suggests a prescription for deriving the ABJM scaling 
function from that of J\f = 4, no such a connection is known for the first and the second 
order coefficients and f^ 2 - 

From the structure of the BDS-like ansatz it is easy to realize that as long as we are 
interested in the non-trivial part of the amplitude we can forget about f^) 2 that would 
contribute only to constants. The lack of information about /<ig\, instead, leaves the 
1/e pole undetermined 2 . Therefore, the two relations above are sufficient to formulate 
an almost complete prediction, up to scheme-dependent subdivergent terms. 

Combining the two ingredients (5.1, 5.3), we can write 

M<£ (^a) = \ (M2(2e,//) + Die))" -6( 2 (^«(4e,/) + D{2e)) 

-(s^ + t^)^ + C^(e) (5.4) 

where C^(e) includes finite contributions coming from and f^ 2 . 

Alternatively, we can collect the pieces that reproduce the Af = 4 SYM amplitude 
at two loops [19] and cast the previous expression into the form 

Mf} (e,/i A ) = .MgW") +M^(2e^")D(e) - 5 C2 M$ (4e, //) 

1 f (4) 
-6C 2 ^(2e) + l -D{ef - (s^ + t^) ^ + C&\e) 

(5.5) 

where fcsi anc ^ ^cs nave been defined so as to include 0(e^ 1 ) and 0(e°) terms re- 
spectively, arising when reconstructing the four dimensional four loop amplitude 3 . 



2 To restrain the lack of information coming from the unknown fcsi, one could switch to a scheme 
where the choice of the regularization scale would be analogous to that in (2.2) (p? A — V \J\ = 2 2e fj,\), 
allowing for no 1/e poles in the series expansion of the two-loop amplitude. In that case fc S \ would 
only affect the coefficient of the 1/e pole at four loops, but not the finite part. 

^ 3 More precisely jgr = - 2 /^i = - | Cs, and 6^ = C<g - 4 f% fl = 

&CS - I C3- 
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Eq. (5.5) expresses our prediction for the four-loop four-point amplitude in ABJM 
in terms of the two-loop amplitude in J\f = 4 SYM theory. An explicit expression for 
its divergent and finite parts is given in Appendix E. In the C coefficient, eq. (E.2), 
we can recognize terms proportional to L 4 , (log 2 2)L 2 and ( 2 L 2 which were present in 
the finite reminder computed in [4]. 

We stress that the most non-trivial contribution to the amplitude, featuring the 
highest weight harmonic polylogarithms, is enclosed in the four dimensional two-loop 
contribution (the first term on the RHS of this equation). Moreover, all the remain- 
ing polylogarithmic dependence on the kinematic invariants is captured by the four- 
dimensional lower order amplitudes. This pattern will show up at any perturbative 
order. 

6. Discussion 

In this paper we have algebraically derived an exact relation between the two-loop 
four-point amplitude of the ABJM theory divided by its tree level expression and the 
corresponding one-loop ratio in the M = 4 SYM theory. That a relation were to be 
at work was already observed in [2, 3, 4] for the two amplitudes computed up to finite 
terms in the IR regulator. Here, we have extended this relation to an exact identity 
holding at any order in e. 

Besides the technical reasons underlying this identity which can be inferred from our 
proof, it would be very interesting to understand whether it has a deeper explanation 
based on more robust conceptual grounds. Moreover, still in this direction, it would 
be important to investigate whether it is an accident of the four-point amplitude or it 
holds also for higher-point amplitudes. 

Assuming that, as in the M = 4 SYM theory, the ABJM amplitudes have an 
iterative structure encoded in a BDS-like exponentiation ansatz, the all-order-in-e 
result at lowest order in perturbation theory would be sufficient for identifying all 
higher order contributions. In three dimensions this is indeed two loops. According to 
our identity, the two-loop 3d amplitude can be rewritten in terms of the lowest order 
contribution to the four dimensional one which enters the ordinary BDS equation for 
the whole Ad amplitude. Therefore, this allows to conjecture that an exact relation 
should hold between the complete four-point three dimensional amplitude and its four 
dimensional counterpart. 

A strong non-trivial test of this conjecture would come from a direct evaluation 
of the ABJM amplitude at four loops, either by traditional perturbative methods or 
by generalized unitarity cuts. The result could in fact confirm or kill the prediction 
that we have for the amplitude at this order and, at the same time, could give some 
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indication on the explicit expression for the scheme dependent coefficient /^i that our 
ansatz leaves undetermined. 

In Refs. [2, 3, 4] a duality between the two-loop four-point amplitude and a 
3d light-like four-polygon Wilson loop [23, 24] was pointed out, as long as the two 
quantities are evaluated up to finite terms in e. It would be interesting to investigate 
whether this duality persists at sub-leading orders. Given our algebraic identity and 
the fact that the Ad amplitude is known up to order e 4 [10], the explicit result for 
the two-loop ABJM amplitude at that order in e is now available. It would then be 
interesting to try and push the evaluation of the Wilson loop up to the same order and 
compare the two results. 
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A. ABJM theory: Notations and conventions 



In M = 2 superspace, the physical content of U(N)k x U(N)_k ABJM theory [1] is 
organized into two vector multiplets (V, V) in the adjoint representation of the first and 
the second {7(iV)'s respectively, and four chiral multiplets A 1 and Bi, i — 1,2, with A 1 
in the (N, N) and Bi in the (N, N) (anti)bifundamental representations. 
The M = 6 supersymmetric action reads [20, 21] 

S = S CS + Smat (A.l) 

with 

K 



>cs 



47T 



d 6 x<f9 / \Tr 



VD a (e- tv D a e tv ) 



— Tr 



VD a f e~ tv 



D a e 



tv 



+ 



t = y rf^rf 4 ^ Tr (A^A'e-v + B^^B^'^ 

— I d*xd 2 6e, k e 3l Tr(A l B ] A k B l ) + — I d 3 x d 2 6 e ik tji Tr(AiB j A k B l ) 



(A.2) 

where K is an integer, as required by the gauge invariance of the effective action. In 
the perturbative regime we take A = 1. 

We are interested in four-point scattering amplitudes. Without loosing generality 
we consider chiral superamplitudes of the type (A l BjA k B{), as the other superampli- 
tudes can be obtained from these ones by SU (4) R-symmetry transformations. 

At any order in perturbation theory they can be inferred from the corresponding 
contributions to the effective superpotential [3, 4]. When going to components, they 
give rise to amplitudes for two scalars and two fermions. 

We write the amplitude divided by its tree level expression as an expansion in 
powers of the 't Hooft coupling 

M = ^\ l M {l) (A.3) 
i 

where A is the dimensionless effective coupling in d = 3 — 2e, and = A^ loops ^ / A tree 
includes the mass scale of dimensional regularization (/i 2 ) /e . 

B. Functions appearing in the result at order e 2 

Here we review the definition of harmonic polylogarithms H ai ___ an which are ubiquitous 
in solving Feynman integrals. These functions are defined recursively as 



H ai a 2 ...a n (x) = I dt f ai (t)H a2 ... an (t) 




(B.l) 
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where 

f±i(x) = T l—, fo(x) = -, (B.2) 

-L — (— Ob Ob 

H ±1 (x) = T\og(lTx) ) H (x)=\ogx, (B.3) 

and the index n is referred to as the weight of the given harmonic polylogarithm. Up 
to weight 4, these functions (when only or +1 indices are present) may be expressed 
in terms of ordinary polylogarithms. 

All the harmonic series arising from closing our MB integrals can be solved ana- 
lytically. Simple and nested harmonic sums are conventionally defined as [22] 

3=1 J 3=1 J 

Here we list some nested series, which are encountered during the order e 2 evaluation 
of integrals, producing high order harmonic polylogarithms 

00 n 1 

Si(n - If- = -2S 1>2 (z) + H^Ho^z) + -H^zf (B.5) 

n=l 

00 n i 

Si(n - l)S 2 (n - 1)- = --%W 2 - #i(z) (5 1)2 (z) - if ,o,i(*)) + 

n=l 

^{zfH^zf (B.6) 
J2s 2 (n- if- = H mi (z) + 2^,0,1,0,1 (z) (B.7) 

n=l 

These are relevant for summing series of second order in polygamma functions ip^(n) 

^(x) = ^] 0g T(x) (B.8) 
such as J2n=i [^ (0) M] 2 ZT=i ^ (0) M* n and XT=i M] 2 respectively 



C. Explicit Comparison 

In this section we compute explicitly the three dimensional amplitude at order e 2 and 
compare it to the four dimensional one. 

In order to make closer contact with the four dimensional case we choose the IR 
scale to be fi" 2 = 47re~ 7E /i 2 . This choice has the advantage that many non-trivial 
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terms multiplied by powers of log 2 drop off the series. The price to pay is to reinstate 
a e _1 pole, which could have been removed from the coefficients Aj by using the scale 
// in (2.2). We write 



J A a 



^(^") = (^ + ^)£f (C.l) 

J=-2 

By e-expanding the integrals arising in the computation and using the tools reviewed 
above we determine explicitly the coefficients. With the shorthands x = s/t and 
L = log x, they read 



A 



2 



f 

4 



A 1 = ~\og2 

1 7T 2 

A = +-L 2 + y + log 2 2 

185C(3) 4 log 3 2 7 2l n 

+ — vr 2 log 2 

12 3 12 s 

A_ 2 = +#1,0,0,1 ("i) + #0,0,1,1 (-^) + #0,1,0,1 ("^) +^0,0,0,1 ("£) 

+#o,i,i {-^L + H lfltl 

>B) i3 -^)-^-B)+K 

7T 

24 

91C(3)log2 233vr 4 4 log 4 2 17 2l 2o rt , 

In order to make a comparison, we report the corresponding result for the one- 
loop four-point amplitude in M = 4 SYM as can be found in Ref. [10]. Using the 
results of Appendix B of that reference, suitably adapted to our conventions (euclidean 
signature, x — > 1/x and e — > 2e) and choosing a slightly different definition for the 
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series expansion 



M2M") = (^ £ + r 2 ^| (C.3) 

i=-2 



we find 

c 2 



1 

d = 



c -' = -^(-;)-^B) L 

+ y Cs 

C- 2 = +#i,o,o,i (-i) + #0,0,1,1 (~) + #0,1,0,1 ("^ J + #0,0,0,1 

+ #0,1,1 + ("^ 

1 / 1\ .3 1 2 _. / l\_ . „ / l\ I I 



Comparing them with eq. (C.2) we see that the two sets of coefficients differ by con- 
stants and by a non-trivial term — ^7r 2 L 2 = — |£ 2 log 2 (s/t) at order e 2 . Owing to this 
observation we conclude that up to order e 2 the following relation holds 

= M$(2e) + £>(e) + O(e') (C.5) 

where A(e) = 1 - 5C 2 e 2 + 0(e 3 ) and D(e) is given in (D.4). This precisely meets the 
expectations from (3.27), which proves this formula to be valid to all orders in e. 

D. Series expansions 

Here we list the e-expansions of the functions used in the computations of Sections 3.2, 
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4 and 5 



A(e) 



r(l - 2e)r(l - 3e)r(l - 4e) 
r 3 (l - e)r(l - 6e) 

l-5C 2 e 2 -40C 3 e 3 - 



2 ,n,J 821 il e 4 + (9(e 5 ) 



(D.l) 



B{e) = 



r 2 (-2e)r(i + 26) 27E£ 

r 3 (l-4e) 

_L_i C2 _| C3£ -| C4£W) 



(D.2) 



+ r-6C 2 log2 + ^C3-^log 3 2)e 



+ (-12 C 2 log 2 2 - ^C4 + y Ca log2 + 1 log 4 2 ) e 2 + O (e 3 ) 



(D.3) 



£>(e) = (*"* + *-*) |-llog2+ (log 2 2-^C 2 ) 
+ (-6C 2 log2-ic 3 -^log 3 2V 



+ (-12C 2 log 2 2-^ 



C 4 + f C 3 log2 + llog 4 2) e 2 + 0( e 3 ) 



(D.4) 



E. The ABJM four— point amplitude at four loops 

In this Appendix we give the explicit result for the ABJM four-point amplitude at four 
loops as derived from eq. (5.5). Writing 

A4£>(e, = + + G(e) (E.l) 

3=0 
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we find 

C 4 



16 

C 3 = ^log2 
C 2 = -^-ilog 2 2 
C l = - A H^ l {-x)-- A H Qtl {-x)L 



+ l H l (-x)L 2 + \it 2 H l (-x) + (x<* -) - I log2L 5 

o o \ X J Z 

65 A 1 1 3 o , 15 A 1 o ^os!l 

~ TT7 Cs - o lo S 2 + — C 2 log2 - 



24^° 3 ° 4 ~ ° 16 



1111 

Co = - 2 #0,0,0,1 (-^ - 2 ^0,0,1,1 (-^ - 2 #i,o,o,i(-z) - ^ #o,i,o,i (-&) 

+ I # ,i,i(-*) ^ + ^ ^iAi(-^) ^ - ^ ^ 
+ I ^ ,i(-x) - | C 2 ^i,i(-x) + \ C 3 ^i(-x) 
+ ^ log 2 (6 tfo,o,i (-*)- 6L#o,i (~x) 

+ 3H 1 (-x)L 2 + 3n 2 H 1 (-x)^j + (^x <* 

+ ll L 4 + 3 2 2 1 2 

24 4 s 4 s2 

+ ±§ C 4 - ^ Ca log 2 + ^ G log 2 2 + | log 4 2 

-ilog2 / S 1 -l / S 2 + ^4 4 s ) (0) (E.2) 
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